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Abstract 

We report computational results of blood flow through a model of the human aortic arch and 
a vessel of actual diameter and length. On the top of the aortic arch the branching of the arter- 
ies are included: the subclavian and jugular. A realistic pulsatile flow is used in all simulations. 
Calculations for bifurcation type vessels are also carried out and presented. Different mathemat- 
ical methods for numerical solution of the fluid dynamics equations have been considered. The 
non-Newtonian behaviour of the human blood is investigated together with turbulence effects. 
A detailed time-dependent mathematical convergence test has been carried out. The results of 
computer simulations of the blood flow in vessels of three different geometries are presented: for 
pressure, strain rate and velocity component distributions we found significant disagreements be- 
tween our results obtained with realistic non-Newtonian treatment of human blood and the widely 
used method in the literature: a simple Newtonian approximation. A significant increase of the 
strain rate and, as a result, a wall shear stress distribution, is found in the region of the aortic 
arch. Turbulent effects are found to be important, particularly in the case of bifurcation vessels. 

PACS numbers: 47.50. Cd, 47.63.Cb 

Keywords: Fluid Dynamics, Navier-Stokes Equation, Human Blood, Aortic Arch, Pressure, Strain Rate, 
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I. INTRODUCTION 



Cardiovascular diseases, such as ischemic heart disease, myocardial infarction, and stroke 
are leading causes of death in Western countries. All of these vascular diseases share a 
common element: atherosclerosis. They also share a common final event: the failure or 

nn 

destruction of the vascular wall structure l2j. 

Atherosclerosis reduces arterial lumen size through plaque formation and arterial wall 
thickening and it occurs at specific arterial sites. This phenomenon is related to hemody- 
namics and to wall shear stress (WSS) distribution [3i]. From the physical point of view WSS 
is the tangential drag force produced by moving blood, and it is a mathematical function of 
the velocity gradient of blood near the endothelial surface ^\. Arterial wall remodeling is 
regulated by WSS j^. For example, in response to high shear stress arteries enlarge. 

Further, from the biomechanical point of view one can conclude, that the atherosclerotic 
plaques localize preferentially in the regions of low shear stresses, but not in regions of 
higher shear stresses. Furthermore, decreased shear stress induces intimal thickening in 
vessels which have adapted to high flow. 

Final vascular events that induce fatal outcomes, such as acute coronary syndrome, are 
triggered by the sudden mechanical disruption of an arterial wall. Thus, we can conclude, 
that the final consequences of tragic fatal vascular diseases are strongly connected to me- 
chanical events that occur on the vascular wall, and these, in turn, which are likely to be 
heavily influenced by alterations in blood flow and the characteristics of the blood itself. 

Currently researchers in the field of biomechanics and biomedicine conduct laboratory 
investigations of human blood flow in different shape and size tubes, which are designed 
to be approximate models of human vessels and arteries, see for example 0, 0]. Some 



researchers also carry out intensive computer simu 
see for example 
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ations of these bio-mechanical systems. 
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Also, in some special laboratory works specific stents are incorporated in such artifical 
vessels (tubes). Stent implantation has been used to open diseased coronary blood vessels, 
allowing improved perfusion of the cardiac muscle. Used in combination with drug therapy, 
vascular repair and dilation techniques (angioplasty) the use of metallic stents has created a 
multibillion dollar industry. Stents are commonly used in many different blood vessels, but 
the primary site of deployment is in diseased coronary arteries. 



Stents represent a very special case in the modeling research problems mentioned above 
M.Ta.i„gi„.oaccou,. .hat ste„ts have a very small size and rather complicated structure 
and shape, this situation makes it difficult to obtain precise measurements. Therefore high 
quality and precise computer simulations of blood flow through vessels with implanted stents 



would be most useful [191]. Work o 



mention several pertinent studies 
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this type is already underway and we would like to 
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Nevertheless, there are still many difficulties in obtaining precise realistic geometries 
for the required vessels. Human arteries, especially the aorta, have complicated spatial- 
geometric and characteristic configurations. For example, the aortic arch centerline does 
not lie on a plane and there are major branches at the top of the arch feeding the carotid 
arterial circulation. One of the main problems in the field of bio-medical blood flow simula- 
tion is to obtain precise geometrical-mathematical representations of different vessels. This 
information in turn needs to be included in the simulation programs. 

However, in our opinion, as a first step of these investigations it would be useful to apply 
simpler 3D-geometry forms and models, but to take into account as much as possible the 
precise physical effects of blood movement such as the non-Newtonian characteristics of 
human blood, realistic pulsatile flow, and possible turbulent effects. Because of pulsatility 
the last one may be significant in contributing to the final results of this study. 

In the current work we carried out real-time full- dimensional computer simulations of a 
realistic pulsatile human blood flow in actual size vessels, vessels with a bifurcation, and in a 
model of the aortic arch. We take into account different physical effects, such as turbulence 
and the non- Newtonian nature of human blood. The next section presents the mathematical 
methodology and the physical model used in this work. The general purpose commercial 
computational fluid dynamics program FL0W3D is used for its basic functionality, but we 
supplemented its capability by adding our routines to obtain the results presented in this 
work. 

Sec. 3 presents results for three vessels of different geometries. The CGS unit system is 
used in all simulations, as well as for presentation of the results. Conclusions and discussion 
comparing our results to well respected studies is included in Sec. 4. 
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II. MATHEMATICAL METHODOLOGY AND PHYSICAL MODELS 



As wc mentioned above, we undertook pulsatile human blood flow simulation experiments 
using different size and shape human vessel arteries. For each spatial configuration one needs 
to provide a specific approach for the numerical solution to the complicated second order 
partial differential equations of the fluid dynamics (FD). For example, for simple cylindrical 
vessels we used the cyhndrical coordinate system: r — (r, ^,0). However, for the aortic 
arch or bifurcated vessels, where there is no cylindrical symmetry, we applied the Cartezian 
coordinate system: r = {x,y,z). In the cases of the aortic arch and bifurcated vessels we 
used up to five blocks of matched Cartezian coordinate subsystems. 

Below we present the FD equations in a general form, because, for each of the special 
cases, considered in this work and the chosen coordinate system, the partial differential 
equations of the fluid dynamics may look different, although we understand, that the general 
differential operator form of the equations is always unique. 



A. Equations 

The equations of motion for the fluid velocity components {u, v,w) in the 3-coordinate sys- 
tem are the Navier-Stokes equations with specific additional terms included in the FL0W3D 
program: 

du 1 . . du . du . du. ^AyV^ 1 dp ^ 
ot Vf ox ay az xVf p ox 



pVf 
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dv 1 . . dv . dv ^ dv. ^A^uv Rdp ^ 

Hi + ("•^•sJ + "^"^ + «^ = -7% ^ + ^' " 

^^{v - - S ■ Vs) (2) 
dw 1 , , dw , -r^dw , dw , Idp ^ ^ , 

IF + T^f"^' * + "^'^ + "'^'^^ = -p & + + ^' - - 

^^{w -w,tj-S -Ws). (3) 

Here, {u,v,w) are the velocity components in coordinate directions {x,y,z) respectively. 
For example, when Cartesian coordinates are used, R—1 and ^ = 0, see FL0W3D manual 



2J]. is the fractional area open to flow in the x direction, analogously for Ay and A^. 
Next, Vf is the fractional volume open to flow, R and C, are coefficients which depend on 
the coordinate system: {x,y,z) or {r,9,z), p is the fluid den sity , Rsor is a mass source 



term. Finally, {Gx,Gy,Gz) are so called body accelerations |2J], {fx,fy,fz) are viscous 
accelerations, {bx,by,bz) are the flow losses in porous media or across porous baffie plates, 
and the final term accounts for the injection of mass at a source represented by a geometry 
component. Mathematical expressions for the viscous accelerations {fx, fy, fz) are presented 
in the Appendix. 

The term = (u^,f^„,iy^) in equation ([1]) is the velocity of the source component. 



which will generally be non-zero for a mass source at a General Moving Object (GMO) |24|. 
The term Us = {ug, Vs, Ws) is the velocity of the fiuid at the surface of the source relative to 
the source itself. It is computed in each control volume as 

'~Ps dA 

where dQ is the mass fiow rate, ps fiuid source density, dA the area of the source surface in 
the cell and n the outward normal to the surface. The source is of the stagnation pressure 
type when in equations. ([HIS]) 6 = 0.0. Next, 6 = 1.0 corresponds to the source of the static 
pressure type. 

As we already mentioned, in all simulations we considered the blood fiow as a pulsatile 



fiow. The final result for the infiow waveform has been taken from figure 3 of work [25| . The 
values of the velocity waveform, which were used in our simulations, are shown in Table 1 
depending on pulse time. The pulse was applied for 5.5 cycle times in our work. 

These velocity values are used as time-dependent infiow initial boundary conditions. 
These numbers are included directly in the FL0W3D program. 

Next, the general mass continuity equation, which is solved within the FL0W3D program 
has the following general form: 

Rdif is a turbulent diffusion term, and Rgor is a mass source. The turbulent diffusion term 
is 

R..f = g^iv^Ax^) + Rg^MR-) + g-^iv,A-) + (6) 



where the coefficient Vp = Cpfi/p, fi is dynamic viscosity and Cp is a constant. The Rgor 
term is a density source term that can be used to model mass injections through porous 
obstacle surfaces. 

Compressible flow problems require the solution of the full density transport equation 
In this work we treated blood as an incompressible fluid. For incompressible fluids, p = 
constant and the equation (2) becomes the following: 

^M.) + ^{vA,) + ^{wA^) +^^ = ^, (7) 
ox ay oz x p 

It is assumed, that at a stagnation pressure source fluid enters the domain at zero velocity. 
As a result, pressure should be considered at the source to move the fluid away from the 
source. For example, such sources are designed to model fluid emerging at the end of a 
rocket or the simple deflating process of a balloon. 

In general, stagnation pressure sources apply to cases when the momentum of the emerg- 
ing fluid is created inside the source component, like in a rocket engine. 

At a static pressure source the fluid velocity is computed from the mass flow rate and the 
surface area of the source. In this case, no extra pressure is required to propel the fluid away 
from the source. An example of such a source is fluid emerging from a long straight pipe. 
Note that in this case the fluid momentum is created far from where the source is located. 

Turbulence models can be taken into account in FL0W3D. It allows us to estimate 
the influence of turbulent fluctuations on mean flow quantities. This influence is usually 
expressed by additional diffusion terms in the equations for mean mass, momentum, and 
energy. The turbulence kinetic energy per unit mass, q, is the following 

where P is shear production, G is buoyancy production, Diff is diffusion, and D is a 
coefficient 2^ . 

When the turbulence option is used, the viscosity is a sum of the molecular and turbulent 
values. For non-Newtonian fluids the viscosity can be a function of the strain rate and/or 
temperature. A general expression based on the Carreau model is used in FL0W-3D for 
the strain rate dependent viscosity: 

^ , /ip-E'r - Poo \2 , . 

^ Aoo + [Ao + (AiET)2e,,e,,](i-)/2 ^ y(e,,e,,) ' ^ > 



where Cij = l/2{dui/dxj + duj/dxj) is the fluid strain rate in Cartesian tensor notations, 
A'cx), A*o, Ao, Ai, A2 and n are constants. Also, Et — exp[a{T*/{T — b) — C)], where T*,a,b, 
and c are also parameters of the temperature dependence, and T is fluid temperature. This 
basic formula is used in our simulations for blood flow in vessels and in the aortic arch. For 
a variable dynamic viscosity /i, the viscous accelerations have a special form. It is shown in 
the Appendix. 

The equations of fluid dynamics should be solved together with specific boundary con- 
ditions. The numerical model starts with a computational mesh, or grid. It consists of a 
number of interconnected elements, or 3D-cells. These 3D-cells subdivide the physical space 
into small volumes with several nodes associated with each such volume. The nodes are 
used to store values of the unknown parameters, such as pressure, strain rate, tempera- 
ture, velocity components and etcetera. This procedure provides values for defining the flow 
parameters at discrete locations and allows speciflc boundary conditions to be set up. 

III. NUMERICAL RESULTS 

Results of our simulations are presented below. One of the most important preliminary 
testing tasks is to check for numerical convergence. This test has been successfully accom- 
plished in this work. A portion of the test calculation results are shown bellow in this paper. 
Next, in this work particular attention has been given to the calculations of wall shear 
stress distribution (WSS). WSS is the tangential drag force produced by moving blood. It 
is a mathematical function of the velocity gradient of blood near the endothelial surface: 
Tw — fJ'[dU{t,y,Rv)/dy]y^Q. Here /i is the dynamic viscosity, t is current time, U{t,y,Rv) 
is the flow velocity parallel to the wall, y is the distance to the wall of the vessel, and Ry 
is its radius. It was shown, that the magnitude of WSS is directly proportional to blood 
flow/blood viscosity and inversely proportional to the cube of the radius of the vessel, in 
other words a small change of the radius of a vessel will have a large effect on WSS. 

A. Straight vessel: cylinder 

First, we chose a simple vessel geometry, that is we considered the shape of a straight 
vessel to be a tube. In our simulations involving a straight cylinder type vessel we applied 
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a cylindrical coordinate system: (r, 9, Z) with the axis OZ directed over the tube axis. 
Different quantities of cells have been used to discretize the empty space inside the tube. 
In the open space (inner part of the tube) the fluid dynamics equations have been solved 
using appropriate mathematical boundary conditions. The size of the tube is: L = 8 cm 
(in length) and R = 0.34 cm (length of inner radius). The thickness of the vessel wall is 
s = 0.03 cm. We have applied 5.5 cycles of blood pulse. 

Let us now consider in more detail the expression ([9]). In these calculations we follow the 
work [9], where the Carreau model of the human blood has also been used. In consistence 
with (9| we choose the following coefficients: A2 = Aqo = 0, a = and = 1, that is we 
don't take into account the temperature dependence of the viscosity. This investigation can 
be done in our subsequent work. Next: Aq = 1, Ai = 3.313 sec, /Xoo = 0.0345 P, /xq = 0.56 P, 
and n = 0.3568. 

Time-dependent results for pressure, strain rate and velocity components V and W are 
presented in Fig. 1. The turbulent effects are not taken into account. We chose to present 
only one precise geometrical point for comparison purposes: the middle point: r = 6 = 0, 
and Z = 4.0 cm. The data for Fig. 1. were obtained with the non- Newtonian model of 
human blood. We refer the reader to the comments provided for the figure. We were able to 
closely replicate the values for all previous cell sizes 2J] and obtain almost identical values, 
for example for pressure, wall shear stress and other parameters, for 0.065 and 0.062 cell 
sizes {2^. This means, that the convergence has been achieved. 

Next, it would be very interesting to compare the results calculated with and without 
the turbulent effect. To support this endeavor we use the realistic non-Newtonian model of 
blood viscosity, the pulsatile flow, and the size of computation cells at which convergence 
has been achieved, that is the 0.062 size for all computational cells 24 1. 

The results are presented in Fig. 2. As we see from Fig. 2 the effect of the turbulence 
is significant, particularly in regard to dynamic viscosity and strain rate. This result means 
that in the case of pulsatile fiows and non-Newtonian viscosity the turbulent term should 
be taken into account. In Fig. 3 we separately show the results for the pulsatile pressure 
distribution and the turbulent energy, again using the middle point of the cylinder. Finally, 
it would also be very interesting to make a comparison between the results calculated using 
both a Newtonian and non- Newtonian viscosity. However, as in previous simulations, we will 
apply the pulsatile fiow with the turbulence included, since it has proved to be important. 
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The results are shown in Fig. 4. As one can see, we obtain significant differences between 
these two calculations. We specifically observed that for the pressure distribution, dynamic 
viscosity and turbulent energy, we obtained significant disagreements. 

Thus, we arrive at the important conclusion: within a time- dependent (pulsatile) flow of 
human blood it is necessary to take into account turbulence and non-Newtonian viscosity. 



B. Hemodynamics in the coronary bifurcation 

Below we show the result of our subsequent simulation involving a 90° bifurcated coronary 
artery in Figs. 5 and 6. The geometrical model of the bifurcation consisted of a 90° 
intersection of two cylinders. This model represents the bifurcation between the left anterior 
descending coronary artery and the circumflex coronary artery. In our opinion, in the case of 
pulsatile flow it is more interesting to present results in a time- dependent way. This method 
can provide a wider picture of highly non- stationary flow systems. In this paper, because of 
space limitations, we just included time-dependent results for pressure, dynamic viscosity, 
turbulent energy, and strain rate. However, we understand, that results which depend on 
spatial coordinates (r, 9, Z) for a few fixed moments of time are also highly useful. 

In the case of the bifurcation shown in Fig. 5, we report the results for only two spatial 
points, which are the two outflow sides: the far right side and the farthest upper side of the 
bifurcation. The length of the lower horizontal vessel is 4 cm and its diameter is 0.54 cm. 
The length of the upper vertical vessel is 1.2 cm and its diameter is 0.4 cm. These sizes are 
consistent with average size human vessels. 

In Fig. 5 blood flows in from the left to the right with the imposed initial velocity profile 
taken from Table 1. The pressure, strain rate and turbulent energy distributions are shown 
for only one specific time moment t=4.329 s. The velocity vectors are also shown on these 
plots. 

Further, Fig. 6 represents our time-dependent results for the two outflow sides mentioned 
above. These results are for pressure, dynamic viscosity, turbulent energy and strain rate. 
The bold black lines are the results for the right outflow side, and the red dashed lines are 
the results for the farthest upper side (see comments to Fig. 6). In conclusion, the main goal 



of these calculations is to adopt them to investigate a case in which a stent 
in the bifurcation area. 



19| is implanted 



in 



C. Blood flow in aortic arch 



The geometry of the blood simulations inside the human aortic arch is shown in Fig. 7. 
On the top of the aortic arch three arteries are included. These arteries deliver the blood to 
the carotid artery and then to the brain. This configuration only models and approximately 
represents the real aortic arch. One of the goals of our simulations is to reveal the physics 
of the blood flow dynamics in this important portion of the human cardiovascular system. 

The aortic arch is represented as a curved tube. The outer radius of the tube is 2.6 
cm. A straight vessel (tube) is also merged to the arch. The length of the straight tube is 
about 4 cm. Again, the thickness of the wall is 0.03 cm, and the inner radius of the tube is 
r = 0.34 cm. Once again we are using the Cartesian coordinate system. We also carried out 
a convergence test. To better represent the shape of the arch we applied five Cartesian sub- 
coordinate systems in our FL0W3D simulations. After the discretization the total number 
of all cubic cells reached about 900,000. It is important to mention here, that we again 
obtained full numerical convergence. In this work we computed pressure, velocity and strain 
rate distributions in the arch, while the human blood is treated as a non-Newtonian liquid 
and while the realistic pulsatile blood flow is used. 

In Fig. 7 we present the results of strain rate distributions inside the arch for two specific 
time moments. At the most left point, which is the inlet, we specify the pulsatile velocity 
source as the initial condition, that is the data from table 1 are used. From the general 



theory of fluid mechanics j4| it is possible to determine together with the blood density and 
viscosity, and spatial geometries, the dynamics of the blood according to the Navier-Stokes 
equation and its boundary conditions. Small vectors indicate the blood velocity. As can be 
seen from Fig. 7 blood flows from left to right in direction. However, because of pulsatility 
blood flows in the opposite direction too. 

The values of the strain rate are also shown. These values are strongly oscillating. From 
the plots one can conclude that in the region of the arch the strain rate values become much 
larger than in the region of the straight vessel. This result represents clear evidence that in 
this part of the human vascular system atherosclerotic plaques should localize less than in 
the straight vessels. However, the higher wall shear stress values in the aortic arch could be 
the reason for sudden mechanical disruption of the arterial wall in this part of the human 
vascular system. These results are consistent with laboratory and clinical observations. In 
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Fig. 8 we depict the pressure distribution in the arch. 



IV. CONCLUSION 

In this work we apphed computational fluid dynamics techniques to support pulsatile 
human blood flow simulations through different shape/size vessels and the aortic arch. The 



realistic blood pulse has been adopted and applied from work [25j. The geometrical size of 
the vessels and the aortic arch have been selected to match the average real values. Human 
blood was treated in two different ways: (a) as a Newtonian liquid when the viscosity of 
the blood has a constant value, and (b) as a non-Newtonian liquid with the viscosity value 
represented by the equation (Q. The numerical coefficients in have been taken from 
work 3. 

It is always difficult to obtain a steady-state cycle profile and stable computational results 
at the very beginning of time-dependent simulations. However, after a short stabilization 
period a steady-state cycle profile can be obtained. In our simulations we used up to 5.5 
pulse cycles to reach complete steady state profiles. We obtained valid results for pressure, 
wall shear stress distribution and other physical parameters, such as the three velocity 
components of blood fiow. All of these were shown in Figs. 1 and 2. 

Our simulations showed that the FL0W3D program is capable of providing stable nu- 
merical results for all geometries included in this work. The time-dependent mathematical 
convergence test has been successfully carried out. Particular attention has been paid to this 



aspect of the calcu 



ations. It is a well known fact that fiuid dynamics equations can have 



fl. 



Therefore, numerical convergence has been tested and confirmed in 



unstable solutions 
this work. 

The result of computer simulations of blood fiow in vessels for three different geometries 
have been presented. For pressure, strain rate and velocity component distributions we found 
signiflcant disagreements between our results obtained with the realistic non-Newtonian 
treatment of human blood and the widely used method in literature: a simple Newtonian 
approximation. 

Our results are in good agreement with the conclusions of works HQ , where the 
authors also obtained signiflcant differences between their results calculated with and with- 
out the non-Newtonian effect of blood viscosity. However, the recent work [26|] should be 
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mentioned, in which the authors performed 2-dimensional simulations of human blood flow 
through the carotid artery with and without the non- Newtonian effect of the viscosity. They 
did not flnd any substantial differences in their results. But in the paper I27I], where the 
authors also performed simulations for the carotid artery, only the non-Newtonian viscosity 
was used. 

Next, the influence of a possible turbulent effect has also been investigated in this work. 
It was found that the effect is important. We believe, that the physical reason of this 
phenomena lies in the strong pulsatility of the flow and in the non-Newtonian viscosity of 
the blood. The contribution of the turbulence is most signiflcant in the area of bifurcated 
vessels. 

Finally, a signiflcant increase of the strain rate and, the wall shear stress distribution, 
is found in the region of the aortic arch. This computational result provides additional 
evidence to support recent clinical and laboratory observations that this part of the human 



cardiovascular system is under higher risk of disruption [28|, |29|] . In future works it would 
be interesting to include the elasticity of the walls of the aortic arch (sol and other vessels. 

In conclusion, we would like to speciflcally point out, that the developments in this work 
can be directly applied to even more interesting and very important situations such as when 
a stent is implanted inside a vessel |l9]. In this case, for example, it would be very useful to 
determine blood flow disturbance, the pressure distribution, strain rate and values of other 
physical parameters. The results of this work should allow us to determine the optimal size 
and shape of effective stents. As we mentioned in the Introduction some research groups 
are carrying out laboratory and computer simulations of blood flow through vessels with 
implanted stents 3]. It is very difficult to underestimate the value of these works. 

V. APPENDIX 



For a variable dynamic viscosity fi, the viscous accelerations are 
d d d ^ 

pVpfx = < - [^{AxTxx) + R^{AyT^y) + Y^{A,T^Z) + - AyTyy)] (lO) 
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(13) 
(14) 
(15) 
(16) 
(17) 
(18) 



In the above equations (fT0|) - (fT2|) the terms w^,Wy and wl are wall shear stresses. If these 
terms are equal to zero, there is no wall shear stress. This is because the remaining terms 
contain the fractional flow areas {A^, Ay, Az) which vanish at the walls 24l |. 
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TABLE I: Time-dependent velocity profile at the vessel inlet. The velocity vector is directed exactly 
over the axis of the tube. Blood velocities Vi are given in cm/sec versus the corresponding time 
moments ti in sec {i = 1, 2, 30). Results taken from Fig. 2 of work 25l |. 
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Figure Captions 



FIG. 1. Test of numerical convergence. Time- dependent dynamic viscosity, strain rate 
and velocity components V and W. Results for a vessel of simple geometry - cylinder type, 
for a specific spatial point inside the cylinder - the middle point. No turbulence effects are 
involved in these simulations with the realistic non-Newtonian viscosity of human blood. 
Black dashed line: calculations with 0.08 size for all cells |24|, red dot-dashed line with 
0.07, green double dot - dashed line with 0.065, and blue bold line calculations with 0.062 
size for all cells. 

FIG. 2. Time- dependent results for a specific geometrical point inside the cylinder: the 
middle point. Black dashed line: simulations without taking into account the turbulence; 
red bold line results with the turbulence. The non-Newtonian viscosity is taken into account. 

FIG. 3. Time-dependent results for pressure and the turbulent energy in the middle 
point of the cylinder. The non-Newtonian viscosity. 

FIG. 4. Results for pressure, dynamic viscosity, turbulent energy and velocity W. 
Time-dependent results for the middle point of the cylinder. Bold black line calculations 
with non-Newtonian viscosity of the human blood; red dashed line with its Newtonian 
approximation. 

FIG. 5. Time-dependent results for a vessel with bifurcation. Pulsatile blood flow, non- 
Newtonian viscosity, and the turbulence effect is included. Bold black line: results for the 
far right outflow side z = 0.0; red dashed line results for the farthest up outflow side y = 0.0. 

FIG. 6. The figures are 2D-plots showing the blood flow in the bifurcated vessels for 
only one precise moment of the discretized time ti = 4.329 sec, the corresponding index is 
i = 40. Upper plot represents the result for the pressure distribution in the bifurcation, 
and the pressure ranges from 2068 dynes/sq-cm to 6758 dynes/sq-cm. The middle plot 
represents the results for the strain rate distribution and the lower plot shows results for 
the turbulent energy in the bifurcation. The range of the values is also shown. 

IS 



FIG. 7. Blood flow in the aortic arch. These two plots represent the full 2D-picture of 
the geometry used in these simulations. Shaded results for the strain rate are also shown, 
the bars on the right show the values. Results are for two specific moments of the time t^Q 
— 4.329 sec and ^41 — 4.440 sec. The values of the strain rate distribution range from 0.0 

1/sec to 357.0 1/sec (upper plot) and from 0.0 to 671 1/scc (lower plot). The maximum 
values of the strain rate are localized in the region inside the arch. Blood flows from right 
to left in both pictures. 

FIG. 8. These two plots represent in more detail the region of the arch together with 
shaded results for the pressure distribution. The bars on the right show the values. These 
results are for two specific moments of the time = 4.329 sec and t4i = 4.440 sec, where 
the pressure ranges from 957 dynes/sp-cm to 5889 dynes/sq-cm (upper plot), and from -616 
dynes/sq-cm to 7070 dynes/sq-cm (lower plot). 
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FIG. 1: Test of numerical convergence. Time-dependent dynamic viscosity, strain rate and velocity 
components V and W. Results for a vessel of simple geometry - cylinder type, for a specific spatial 
point inside the cylinder - the middle point. No turbulence effects are involved in these simulations 
with the realistic non-Newtonian viscosity of human blood. Black dashed line: calculations with 
0.08 size for all cells 2j], red dot-dashed line with 0.07, green double dot - dashed line with 0.065, 



and blue bold line calculations with 0.062 size for all cells. 
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FIG. 2: Time-dependent results for a specific geometrical point inside the cylinder: the middle 
point. Black dashed line: simulations without taking into account the turbulence; red bold line 
results with the turbulence. The non-Newtonian viscosity is taken into account. 
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FIG. 3: Time-dependent results for pressure and the turbulent energy in the middle point of the 
cylinder. The non-Newtonian viscosity. 
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FIG. 4: Results for pressure, dynamic viscosity, turbulent energy and velocity W. Time-dependent 
results for the middle point of the cylinder. Bold black line calculations with non-Newtonian 
viscosity of the human blood; red dashed line with its Newtonian approximation. 
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FIG. 5: Time-dependent results for a vessel with bifurcation. Pulsatile blood flow, non-Newtonian 
viscosity, and the turbulence effect is included. Bold black line: results for the far right outflow 
side z = 0.0; red dashed line results for the farthest up outflow side y = 0.0 
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FIG. 6: The figures are 2D-plots showing the blood flow in the bifurcated vessels for only one 
precise moment of the discretized time U = 4.329 sec, the corresponding index is i = 40. Upper 
plot represents the result for the pressure distribution in the bifurcation, and the pressure ranges 
from 2068 dynes/sq-cm to 6758 dynes/sq-cm. The middle plot represents the results for the strain 
rate distribution and the lower plot shows results for the turbulent energy in the bifurcation. The 
range of the values is also shown. 
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FIG. 7: Blood flow in the aortic arch. These two plots represent the full 2D-picture of the geometry 
used in these simulations. Shaded results for the strain rate are also shown, the bars on the right 
show the values. Results are for two specific moments of the time t^o = 4.329 sec and t^i = 4.440 
sec. The values of the strain rate distribution range from 0.0 1/sec to 357.0 1/sec (upper plot) 
and from 0.0 to 671 1/sec (lower plot). The maximum values of the strain rate are localized in the 
region inside the arch. Blood flows from right to left in both pictures. 
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FIG. 8: These two plots represent in more detail the region of the arch together with shaded results 
for the pressure distribution. The bars on the right show the values. These results are for two 
specific moments of the time t4o = 4.329 sec and t4i = 4.440 sec, where the pressure ranges from 957 
dynes/sp-cm to 5889 dynes/sq-cm (upper plot), and from -616 dynes/sq-cm to 7070 dynes/sq-cm 
(lower plot). 
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